setwd(getwd())

source("libraries.R")
source("study1.R")
source("study2.R")
source("study3.R")
source("study4.R")
source("study5.R")
source("study6.R")
source("study7.R")
source("study8.R")
source("construct_fig2_3_4_S2.R")


## Note: to reproduce Fig S7, Fig S8, Fig S9, Fig S10, Fig S11, Fig S12, Fig S13, Fig S14, and Fig S15 (i.e., the main text figures after excluding outliers)
## you must enter the separate code files and set "exclude_obs <- TRUE" before running the code chunks. 